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SUMMARY 


A necessary condition for neutral oscillatory stability is that 
Routh’s discriminant R, formed from the coefficients of the stability 

equation, is equal to zero. The expression for R iB D(BC — AD) — B 2 E 
where A, B, C, D, and E are the coefficients of the lateral-stability 
equation. In a large number of the cases considered in this stiidy, it 
has been found that the term. B 2 E may be neglected. Routh’s discrimi- 
nant is then factorable into two simplified expressions, that 
is, BC — AD = R^ hnd Dj and either R^ = 0 or D = 0, or both, 

may constitute a condition of neutral stability. Test functions have 
been derived which, if satisfied, indicate that the simplified expressions 
may be used to approximate R = 0. If R 1 = 0 and D = 0 satisfy the 

necessary and sufficient conditions for a neutral-oscillatory-stability 
boundary, D = 0 represents the boundary for the oscillation which has 
a period comparatively longer than the period of oscillation for 
which Rq_ = 0 is the boundary. 

In general, the results of the computations obtained from R^ = 0 

and D = 0 show very good agreement with the results calculated by the 
exact expression for R = 0. The nature of the modes of motion as a 
function of the directional— stability derivative and the effective- 
dihedral derivative is discussed in detail. 


INTRODUCTION 


The results of recent investigations (references 1 and 2 and 
unpublished results of lateral-stability analyses for several experimental 
high-speed airplanes) have indicated that smal 1 variations in some of the 
airplane mass and aerodynamic parameters may cause a pronounced change in 
the oscillatory stability of the airplane. It has been difficult to 
explain the reasons for such pronounced changes because of the complexity 
of the expression for neutral oscillatory stability. This expression, 
based on the lateral— stability equations with three degrees of freedom, 
involves a large number of combinations of the mass and aerodynamic 
parameters. In order to predict the stability of the lateral oscillation. 
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therefore, it appears necessary to make a separate stability analysis 
for each airplane. 

The simplified expressions derived for the neutral— oscillatory- 
stability boundary in the present theoretical investigation simplify 
the calculations required to obtain the boundary in the analysis 
essential for each airplane. Because of the comparative simplicity of 
these expressions, an insight into the important combinations of mass 
and aerodynamic parameters that affect the lateral oscillatory stability 
is also provided. Through further investigation and analysis of the 
effects of these major parameters, the necessity of making separate 
calculations for each airplane might possibly be eliminated. Test 
functions are given which, if satisfied, indicate that the simplified 
expressions may be used. 

The nature of the modes of motion as a function of Cn and , 

P P 

the directional— stability derivative and effective— dihedral derivative, 
respectively, are shown to depend upon the location of the stability 
boundaries plotted as a function of (1, and C 7 . 

0 0 

The results of the calculations based on the simplified expressions 
are presented for comparison with the results obtained by the complete 
expression for the neutral— oscillatory— stability boundary. 


SYMBOLS AHD CGEEFICIEriTS 

\ J I angle of bank, radians 

i]r angle of azimuth, radians 

0 angle of sideslip, radians (v/7) 

v sideslip velocity along the Y-axiB, feet per second 

V airspeed, feet per second 

p mass density of air, slugs per cubic foot 

q dynamic pressure, pounds per square foot (^P^) 

b wing span, feet 

S wing area, square feet 

¥ weight of airplane, pounds 
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m 

g 

T] 


e 


e 
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E 


“X 


% 

*XZ 


mass of airplane, slugs (W/g) 

acceleration due to gravity, feet per second per second 

relative-density factor (m/pSb) 

inclination of principal longitudinal axis of airplane with 
respect to flight path, positive when principal axis is 
above flight path at the nose, degrees (see fig. 1) 

angle between reference axis and horizontal axis, positive 
when reference axis is above horizontal axis, degrees 
(see fig. 1) 

angle between reference axis and principal axis, positive 
when reference axis is above principal axis, degrees 
(see fig. 1) 

angle of flight path to horizontal axis, positive in a 
climb, degrees (see fig. 1) 

radius of gyration in roll about principal longitudinal axis, 
feet 

radius of gyration in yaw about principal vertical axis, feet 


nondimenSional radius of gyration in roll about principal 
longitudinal axis ^kj^/bj 

nondimensional radius of gyration in yaw about principal 
vertical axis ^k^/bj 

nondimensional radius of gyration in roll about longitudinal 
stability axis cos ^ti + E^ 2 sl^T] ^ 

nondimensional radius of gyration in yaw about vertical 


stability axis (VV COS^T| + Ej^ 2 Bin^Tj 


) 


nondimensional product— of— inertia parameter 


(fo 0 2 -%o 2 ) 8in r\ cob T\^j 
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trim lift coefficient 


/W cos y\ 

) 


/Boiling moment \ 
rolling-moment coefficient I J 

yawing-moment coefficient 


lateral— force coefficient 


/Lateral force 


V 


qS 


) 


effective— dihedral derivative, rate of c hang e of rolling- 
moment coefficient with angle of sideslip, per radian 

(acz/ap) 

directional— stability derivative, rate of change of yawing- 
moment coefficient with angle of sideslip, per radian 

(ac n /a P ) 

lateral— force derivative, rate of change of lateral— force 
coefficient with angle of sideslip, per radian 

damping-in-yaw derivative, rate of change of yawing-moment 

coefficient with yawing-angular^-velocity factor, per radian 

^rb 




( 


;e or cnange ox yawmg-momau o uuuiiiuj-tuiu 
angular— velocity factor, per radian ^5C n /c 


rate of change of yawing-moment coefficient with rolling— 

^pb' 

27 , 


damping— in-roll derivative, rate of c h a ng e of rolling- 

moment coefficient with rolling-angular-velocity factor, 

.pb' 


per radian 


rate of change of rolling— moment coefficient with yawing— 
angular-velocity factor, per radian 

rate of change of lateral-force coefficient with rolling- 
angular— velocity factor, per radian 

rate of c hang e of lateral— force coefficient with yawing— 
angular-velocity factor, per radian (SOy/S^) 


time, seconds 
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®b 

Db 

R 

X 


P 

T l/2 


nondimensional time parameter ‘baaed on span (Yt/b) 

differential operator ( 4 — ] 

\ ds b ) 


Routt’s discriminant 


complex root of stability equation 


AX a 


BX 


3 + 


CX 4 - DX + E 0 


complex root of stability equation 
AX’^ + BX’ 3 + CX« 2 + DX* + E = 0 

period of oscillation, seconds 


(X = £ ± ico) 


(X* = i* ± ico’) 


time for amplitude of oscillation to change by factor of 2 
(positive value indicates a decrease to half— amplitude, 
negative value indicates an increase to double amplitude) 


A,B,C,D,E coefficients of lateral— stability equation 


EQUATIONS OF MOTION 


The nondimensianal linearized equations of motion, referred to the 
stability axes, used to calculate the spiral— stability and oscillatory- 
stability boundaries for any flight condition, are: 

Rolling 

2 % (’kW + SxzV*) - + t\V 

Yawing 

^ (^zV* + + + i°n r V ■ 

Sideslipping 

%l, (V + V) = C^p + + (0 L tan ?)* 
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X&b A&u 

When (p 0 e is substituted for <j>, ^e for 

for P . in the equations written in determinant form, A 
root of the stability equation 


and P Q e 
must he a 


AX 


^ + BA 3 + CA 2 + DA + E = 0 


(1) 


where 


V 


A = S^ 3 (Sj 2 ^ 2 - Kxz 2 ) 

B = - ^ 2 (^V% + % 2 % + *z\ - <*iz\ 

~ %Z C 2 r ~ %Z C n p ) 


C= ^x\% + + + 2 V*p - 

- KWip ~ %%;% ” 2 %% + *22%% ~ 

- % 2 Cy r Cn p + 

D = “ i Wp “ ^%% + J % C 2 r % + ^%% 


+ 2tJ h C L K XZ Cn p ~ 2 ^TD°1^Z 2 °lp i ~ 2^ tan y 

+ 2 ^h K xz C 2 p c L tan 7 + | c z %% - J- C^C^Cy^ 


- | C 2r Cn p Oy p + i C^C^Cyp 


E = 


2 %%% ~ %%) + | C L tan 7(0^ 
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Multiplying equation (1) "by p, and substituting X = — results 

h) 

in the stability equation 


AX*^ + BX ' 3 + CX ’ 2 + DX* + E = 0 


( 2 ) 


where 



C 

C = — 


D = D 


1 = I^E 

The damping and period of the lateral oscillation in seconds are 
given respectively by the equations 

_ 0.69 

1/2 " V 

P = £^8 ^ 

a>* V 


where £’ and a>* are the real and Imaginary parts of the complex root 
of stability equation ( 2 ). 
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AHALTSIS 


The necessary and sufficient conditions for neutral oscillatory 
stability, as shown in reference 3 , are that the coefficients of the 
stability equation satisfy Routes discr iminant set equal to zero 


B = BCD - AD 2 - B^ = 0 


(3) 


. and that B and D have the same sign. The expression for R = 0 c an 
be derived by assuming that the stability equation has two roots X = ± ico, 
where cn is the angular frequency of the neutrally stable oscillation. 

This assumption is based on the fact that for the condition of neutral 
oscillatory stability the real part of the complex root must be zero. 

If X = ico is substituted in the stability equation, the following two 
equations are obtained: 



+ E = 0 


W 


— Bcc>3 + Dcd = 0 


(5) 


Solving equation (. 5 ) for co 2 and then substituting the result ( oi 5 = £ 

„ , \ B, 

in equation (t) results in Routh's discriminant 


BCD - AD 2 - B%5 = 0 


It is seen from equation ( 5 ) that a> = 


D , 

— is 


the angular frequency of the 


neutrally stable oscillation only when B and D are of the same sign 
since cd must have a real value if the root X = ± io) is to represent 
an oscillation. If B and D are of opposite sign and R = 0 is 


satisfied, the two roots of the stability equation given by X = ±. ice 

are real, one positive and one negative. It is important to note that 
the A, C, and E coefficients may be of opposite sign to the B and D 
coefficients, and neutral oscillatory stability will still occur as long 
Routh’s discriminant is equal to zero and the D and B coefficients 
are of the same sign. 


In general, the R = 0 boundary in the 
branches. The two branches result from the fact that 


plane has two 
R = 0 is a 
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quadratic equation in C, and thus has two C, roots for every value 

ip ift 


3 


of CL . Usually, the two tranches can. he approximated by simplified 
P 

expressions for B = 0. In certain cases, however, which are discussed 
in the section entitled "Test Functions," either one or none of the 
branches may be approximated. 


Wow, the condition E = 0 is a necessary but insufficient condition 
for neutral oscillatory stability. The simplified expressions, therefore, 
which approximate E = 0 do not necessarily represent boundaries of 
neutral oscillatory stability. Other conditions, elaborated on in the 
section "Validity of D = 0 and E-^ = 0 as Weutral— Oscillatory— Stability 

Boundaries," must be satisfied before either expression represents a 
valid boundary. 


There are, therefore, two kinds of tests that must be made: First, 

tests to determine whether E = 0 may be approximated by simplified 
expressions; and, second, tests to determine which of theBe expressions 
represents a boundary of neutral oscillatory stability. The significance 
of the lateral— stability boundaries is indicated by a discussion of the 
modes of motion in the Cjj ,C^ plane. 

P p 


Derivation of Simplified Expressions 


Two of the most important stability derivatives affecting lateral 
oscillatory stability are the directional— stability derivative Cn and 

• p 


the effective-dihedral derivative (L 


'0 


The boundary for neutral oscillatory 


stability is usually plotted as a function of these two derivatives 
with CL as the ordinate and C 7 as the abscissa. The method used to 

3 P 

obtain the neutral oscillatory stability boundary is first to substitute 
the values of the mass and aerodynamic parameters of a specific airplane 
in the coefficients of the stability equation while maintaining 

0 

and C 7 as variables and then to calculate the C 7 roots of equation ( 3 ) 

ip ip 

for several values of C^. For a given value of C^, the expression 
E = 0 is a quadratic equation in C 7 that is of the form 

p 


a l c Zp 


2 


+ 


Cj = 0 


for 
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For a Ter y large number of cases. It has been found that the last term 
of South’s discriminant B^E, which' contributes only to b-]_ and Cj_, 

produces a negligible change in the expression- 

+ b l C 2p + C X = 0 

\ 

If, therefore, the term B 2 E is neglected, equation ( 3 ) for Routh’s 
discriminant reduces to 


R = D(BC - AD) = 0 ( 6 a) 

Thus BC — AD = 0 and D = 0 are the approximate factors of the 
quadratic equation 


a l C 2 c 


+ ^1 C 2 { 


= 0 


The expression BC — AD is henceforth called R^ 

In order to simplify the expressions for R^ and D, the expected 

range of values of mass and aerodynamic parameters for high-speed air- 
craft were substituted in the coefficients of the stability equation to 
determine which terms could be omitted without appreciably affecting 
the values of R-^ and D. The following simplified coefficients were 

obtained: 

A = - %/) 

B - - + + C 'l p % 2 ) 

0 * “ ll|i b K KZ C !p + I Cn r °Jp + % ~ 5 °n p C I r ) 

E - - *b(°l p ~ 20 lfe) c n p 
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E ■ I - w) 

The expressions for E-j_ and D thus hecome 

E i - (Viz - Vi)% + hz(^h. + 0 Ip %z) - % 2 (\ - c lp % 2 ) c np 
-%<)-%% =0 (6b) 

and 

B - *a% - (°j p " aoifeW - 0 ■ (&) 

where 

*1 = E X G n r + ^X^Z °Yp + c l^t 

SB &. i 

' a 3 - %V - %z 2 

The simplified expressions E = 0 and D = 0, as presented, are 

applicable only to conditions of level flight or to conditions of amen 
angles of glide or climb. Simplified expressions for conditions of 
large angles. of glide or climb can be derived by a procedure similar to 
the one presented. 

Test Functions 

The appro xi m a te discriminants B-i = 0 and D = 0 are based on the 

O 

assumption that B E can be neglected when South’s discr imin ant is set 
equal to zero. Thus, the simplified expressions for the neutral— 
oscillatory-stability boundary, E 1 and D, should not be used if 
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p 

including the terms B E causes an appreciable change in the roots 
of R = 0. In appendix A test functions are derived which indicate the 
incremental change in the roots of Rj_ = 0 and D = 0 due to the 

terms B%. If certain conditions placed upon these test functions are 
satisfied, then R^_ and D can he used to calculate the R = 0 

boundary. 


If, at a given value of , the root of R-, = 0 is denoted 

P -L 

by C, = r, the appro ximat e deviation of this root from a root of R = 0 
is given by 


Ar = 


e-L(e - r) 
r ldl(d — r) + e^ 


( 7 ) 


If Ar is small, then R-^ = 0 is a close approximation to one branch 
of R = 0. A suitable criterion for this approximation iB 


|Ar| ^ 


r 

5 


or 

j Ar | = 0.01 

whichever is the larger. 

Similarly if a root of D = 0 is denoted by = d, the approximate 

deviation of this root from a root of R = 0 is given by 


Ad 


62(0 - d) 

r l^l^ r — ^L) + e 


( 8 ) 


If Ad is small, then D = 0 is a close approximation to cne branch of 
of R = 0. A suitable criterion for this approximation is 



< 


d 

5 


or 


Ad 


= 0.01 


whichever is the larger. 
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The expressions for r ±J d-^ e-j_, r, d, and e for use in 

equations (7) and (8) are 


r l 8 ^ i b(- A l K IZ ^2^) 
d l = ^^2 


e l " 

(%) Er o 


r = 


d = 


e = 


(%) 

' P / D=0 


(%) 


B%=0 


where A-^, A^, and A^ are defined in equations (6d) . 


The value- of to he used in theBe test functions should he 

selected from the probable range of C n ' of the airplane for which the 

0 


lateral— stability analysis is to be made. Thus, the approximation 
of Rj_ = 0 and D = 0 to R = 0 is determined in that region of 


the C-p , C •j 

V *0 


plane pertinent to a particular analysis. 


Validity of D = 0 and R^ = 0 as Neutral— 


OBcillatory-Stability Boundaries 


As mentioned previously, for R = 0 to be a boundary for neutral 
oscillatory stability, the coefficients B and D must be of the same 

sign. The three predominant terms of B contain the factors -C n 

“r 



and — , respectively. For positive damping in roll, 

P 
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is negative; and for positive weathercock stability ^C n positive^, 

Cy^ and are negative. Thus, B is positive in the usual case 

where there is weathercock stability and damping in roll. Therefore, B 
must generally be positive in order that R = 0 be a neutral— stability 
boundary. If the exact boundary R = 0 has been calculated, it is 
merely necessary to plot B = 0 and R = 0 and note whether R = 0 is 
located on the side of B = 0 where B is positive. A primary purpose 
of the present paper, however, is to obviate calculation of the exact 
boundary by the US9 of simplified boundaries together with test functions. 
A method to determine the sign of B from the results of the simplified 
expressions is therefore presented in the following paragraph. 


C n (selected from the probable range of 

P 


For a given value of . _ w . 

P V V 

of the airplane for which the lateral-stability analysis is to be madej 


let d be a value of C- 


'P 


different value for which R = R-jB — B^E = 0. 
of C 7 = d* gives 


for which B = 0 and d* be a slightly 

The substitution 


B(d‘) 


B 2 E(d«) 

R^d*) 


(S = a ’) 


The sign of B at the R = 0 boundary 
determined by the signs of E 
little from d’ , the signs of 
the signs of E(d’) and R-j^d* ), respectively, 
at R = 0 is the same as the sign of E/R^ at 

is. 


is therefore 


and Ry at d* . But since d differs 
E(d) and Ry(d) will be the same as 
Hence the sign of B 
B = 0 (fig. 2); that 


B(d’) 


B(d) 

R x (d) 


If the signs of E and R^ are the same, B is positive and represents 
a neutral-oscillatory— stability boundary; if E and R^ are of opposite 

sign, B is negative and then represents a boundary for which the roots 
of the stability equation are equal and opposite in sign. 

The preced ing analysis is applicable for the large majority of cases 
where E(d’) and R^d 1 ) are of the, same sign as E(d) and R^d), 
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respectively. For these cases, the 
the E and R n curves. If the D 


D curve is widely separated from, 
curve is close to either the E 


or Ri curve, the signs of E and R]_ should he determined at = d* . 

However, a very good approximation to the value of d’ can he obtained hy 

adding to d the value of Ad calculated in the previous section entitled 

”Test Functions.” Hence, the sign of D is determined from the signs 

of E and En at C 7 = d + Ad. 

. x 


If the value of C 


■, at which E 7 = 0 is substituted in D and 
ip x 


the resultant sign is positive, Ej_ = 0 is a neutral— oscillatory- 
stability boundary. 


It is interesting to note that for some aircraft, the D = 0 curve, 
which approximates one branch of the E = 0 curve, is a neutral- 
oscillatory— stability boundary over one section of the curve and a 
boundary for equal and opposite real roots over the remaining section. 
This division of the D = 0 curve into two distinct parts is caused by 
a change in sign of v the D coefficient at some point on the curve. If, 
as has been found in a large number of cases, E^ in positive for al l 

values of CL, and C 7 on the D = 0 curve, the sign of the D coef— 

Up L ft 


'0 


ficient depends only on the sign of E at these points. As shown in 
figure 3, therefore, the point of intersection of the curves D = 0 
and E = 0 is the point of separation of the D = 0 curve into two 
characteristically different sections. For points on the hatched side 
of E = 0, the E coefficient is negative and, therefore, the daBhed 
part of D = 0 is a boundary of equal and opposite real roots. Conversely, 
on the unhetched side of E = 0, the E coefficient is positive and the 
solid part of D = 0 approximates a boundary of neutral oscillatory 
stability. 

For Sinai 1 positive or negative values of and negative damping 

in roll, it is possible for B to be negative. A similar analysis is 
applicable to this case where now D must be negative to satisfy the 
necessary condition that E = 0 is a boundary of neutral oscillatory 
stability. 

In general, when the simplified expressions are used to obtain a neutral- 
oscillatory— stability boundary, the procedure to be used is as follows: 

(l) For a given value of Cm , selected from the probable range 

P 

of C,. of the airplane for which the lateral— stability analysis is to 


x 0 

be made, calculate r and d, the C 7 roots of 

4 0 


E 1 = 0 


and D = 0, 


respectively. 
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(2) Determine the value of Ar and Ad "by substituting the results 
into the test functions. 

( 3 ) If the criterions for Ar and Ad as set forth in appendix A 
are satisfied, consider R}n.= 0 and D = 0 close approximations to 

the R = 0 boundary. 


(*0 In order to determine the validity of R-^ = 0 as a boundary of 
neutral oscillatory stability, substitute the given value of 


and C-; 




= r 


into the D coefficient. (If the resulting sign 1 b 


positive, R^ = 0 approximates a branch of the neutral-oscillatory— 
stability boundary. ) 


(5) In order to determine the validity of D = 0 as a boundary of 
neutral oscillatory stability, substitute the given value of C n 

1? P 

and C ^ q = d into — . (If the resulting sign is positive, D = 0 
P R^ 


approximates a branch of the neutral— oscillatory— stability boundary; if 
the resulting sign is negative, D = 0 approximates a boundary of equal 
and opposite real roots.) 


Nature of Modes of Motion in the C n ,C 7 Plane 

n P 4 P 


In this section^ the changes in the roots of the lateral— stability 
equation, which' occur upon crossing the various stability boundaries, are 
discussed according to the principles of the theory of equations as 
given in references 3 and 4. The solution of the lateral-stability 
equation gives four roots which may be four real roots, two pairB of 
conjugate complex roots, or two real rootB and one conjugate complex pair. 
A pair of complex roots indicates an oscillatory mode and a real root 
indicates an aperiodic mode. If the airplane is disturbed from its 
trimmed condition by an arbitrary disturbance, the subsequent motion is 
compounded of these modes in different proportions. The method of calcu- 
lating the different proportions of the modes is presented, for example, 
in references 5 and 6. Such calculations of the motion for numerous 

points throughout the Cn ,C 7 plane would be very laborious. It is 

P P 

more practical, therefore, to investigate merely the types of modes that 

may be expected throughout the On ,C- 7 0 plane as indicated by the 

P P 

stability boundaries. The calculation of the motion could then be limited 
to several points of interest. 


Consider a case where the neutral— oscillatory— stability boundary R^ = 0 
and the spiral— stability boundary E = 0 are located in the first quadrant of 
figure lj-(a). The area between the two boundaries is a region of complete 
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stability. 

and C, „ 

Z 0 


The roots of the stability equation for combinations Of C4 
in this region, such as point A in figure 4(a), are two negative 


n 0 


real roots and one conjugate complex pair with the real part negative. One 
of the real roots wliieh is numerically small corresponds to the spirally 
stable motion of the airplane. The other real root corresponds to the 
heavy damping of the pure rolling motion. The complex roots with the real 
pert negative show that the so-called Dutch roll oscillation is stable. 
Passing through the E = 0 boundary from point A to point B causes the 
spiral mode to become unstable, and crossing through the B-^ = 0 boundary 

from point A to point C causes the oscillatory mode to become unstable. 

The second branch of the B = 6 boundary plotted in the second quadrant 
as D = 0 is not a neutral-oscillatory-stability boundary but rather a 
boundary for equal and opposite roots as determined by the analysis 
presented in the section entitled "Validity of D = 0 and Bj_ = 0 as 


ITeutral— Oscillatory-Stability Boundaries . " 
equation for combinations of C n ^ and C 2 ^ 


The roots of the stability 
on this boundary are two equal 


and opposite real roots and a pair of complex rootB with the real part 
negative. The positive real root is the spirally -unstable mode, and the 
negative real root is the damping— in— roll mode. The oscillation continues 
to remain stable even though the D coefficient is negative. 


\ 


For the case where one oscillatory— stability boundary D = 0 appears 
in the first quadrant and another oscillatory stability boundary B 1 = 0 

is in the second quadrant, (fig. 4(b)), the period of the neutrally stable 
oscillation is much greater on D = 0 than on B^ = 6. This fact can be 

shown to be true by investigating the angular frequency of the neutrally 
stable oscillation for points located on the B 1 = 0 and D = 0 boundaries. 

As shown previously, the angular frequency cd is equal to \Jd/B; and, 
therefore, since the boundary D = 0 approximates one branch of. B = 0, 
the angular frequency for points on that branch is very small. For 
combinations of G n and on Bj_ = 0, the angular frequency is much 

greater. In general, D = 0 is a neutral— oscillatory— stability boundary 
/for a long— period oscillation. The roots at point A of figure 4(b) have 
the same character as the roots at point A of figure 4(a), that is, two 
negative real roots and one pair of conjugate complex roots. At point J3 
the roots of the lateral— stability equation are two pairs of conjugate 
complex roots. It is interesting to note that the boundary for two equal 
roots occurs between point A and point B and can be considered the boundary 
beyond which two pairs of complex roots exist. Beference 4 shows that for 
a quartic equation 


AX 


+ BX 3 + 


CX + DX + E = 0 
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the "boundary for equal roots is obtained "by setting the discr iminant 


- Ij-P 3 - 27Q 2 


equal to zero, where 



and 


Q 


= _ B 2 E + BTO + 


8 ^ E _^2 


2C 3 

27 


Between this "boundary and D = 0, the period of the stable oscillation 
which corresponds to the newly formed pair of complex roots is longer 
than the period of the oscillation which corresponds to the other pair 

of complex roots. As C ? is increased to point C on the unstable side 

3 

of D = 0, the newly formed long— period oscillation is the one that 
becomes unstable whereas the short-period oscillation remains stable. 

At point D the roots consist of a spirally unstable mode, a stable 


mode due to the derivative C 


unstable in passing through 


Ip* 
El = 


and a stable oscillation which becomes 
to point E. 


Figure h(c) represents the case where both = 0 and D = 0 

appear in the first quadrant but only E^ is a neutral— oscillatory— 

stability boundary. The curve D = 0 is the boundary for two equal 
and opposite real roots. Point A once again has two real negative rootB 
and a pair of complex roots with the real part negative. At point B, on 
the unstable side of B-^ = 0, the real part of the complex roots is 


positive and indicates an unstable oscillation, whereas the two real 
roots are still negative. The boundary for C = 0 is between E^ = 0 

and D = 0. Some Investigators of lateral stability have thought that 
a radical change occurs in the roots upon crossing through this boundary. 
The calculations indicate, however, that the roots do not vary appreciably 
upon passing through C = 0. At D = 0, however, there must exist two 
equal and opposite real roots; this condition is possible only if the 
complex roots divide into real roots since the other two real roots are 
negative in sign. The calculation of roots at point C indicate that the 
complex roots had divided into two real positive roots, one of which was 
exactly equal in magnitude to one of the negative roots. Again, the 
boundary for two equal roots, located between C = 0 and point C, would 
determine the combination of 


(L, and C; where the complex roots 

p ^ R 


'3 


divide into tiro real roots. 
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There have "been several cases for which a neutral— oscillatory— 
stability boundary did not exist in the C n ,Cj plane. An analysis of 

these cases indicated that the boundary for equal roots was in the 
oscillatorily stable region and had divided the stable oscillation into 
two subsiding modes. The neutral-oscillatory— stability boundary, therefore, 
would not have any significance. 


RESUITS MB DISCUSSION 


The simplified expressions were used to calculate R 1 = 0 and D = 0, 

and the results are compared with the results of the calculation of R = 0 
based on the complete expression. Not only do R^ = 0 and D = 0 show 

good agreement with R = 0 (figs. 5 to 13) but the comparative simplicity 
of the Rj and D expressions allows identification of the major 

parameters that affect the stability boundaries. 


Effect of 


S - 


on the Branch of R = 0 


Approximated by D = 0 


Reference 2 shows 
is obtained when C n 

P 


that a stabilizing shift in the ,R = 0 boundary 
is increased in a positive direction up to a certain 


value, but further increases in the positive direction cause a destabi- 
lizing shift in R = 0.. The effect of varying C n ^ on the R = 0 curve 


is presented in figure 5 for a model tested in the Langley free— flight 
tunnel. The figure illustrates very good agreement between R = 0 and 
the simplified expressions R^ = 0 and D = 0. The expression for D = 0 


is 


(% - gc Az)S 

S - 2 °I?2 2 


which indicates that for positive CL. when the numerator ie negative in 

P 

sign, the D = 0 boundary is in the second quadrant for negative values 
of Cup — 2CjK^ — Aq and in the first quadrant for positive values of Ag. 

For the cases of negative Ag presented in figure 5, the D = 0 boundary 

would appear in the second quandrant. It can be shown, however, by the methoc 
described in the section "Validity of D = 0 and Rj = 0 as Neutral- 

Oscillatory— Stability Boundaries" that D = 0 in the second quadrant is not £ 
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neutral— oscillatory— stability boundary and bence is not plotted in 
figure 5 . However* as is increased in a positive direction* where 

now A 2 is positive* an increase in the positive value of Ag causes 
the D = 0 boundary to shift upward in the first quadrant in a 
destabilizing direction. 

From, the results shown in figure 5* it is seen that for the cases 
of Cn^ equal to 0.30 and 0.40 only the solid— line part of 

the R = 0 curve in the first quadrant (the branch which may be approxi- 
mated by D = 0) is a neutral— oscillatory— stability boundary. The short- 
dash— line part of R = 0 is a boundary of equal and opposite real roots. 
The reason for this division of the R = 0 curve into two parts is 
discussed in the section entitled "Validity of I) = 0 and R^ = 0 as 

Neutral— Oscillatory-Stability Boundaries" and is illustrated in figure 3 . 


Effect of — 2 CjK 2 on the Branch of R = 0 

Approximated by R^ = 0 


The important effect of C n on R^ = 0 occurs only in the coef- 
ficient of the term* A 1 K 77 — AgA^* in which affects only 

the factor Ag. The sign of is always positive and the sign of A^ 

is negative for positive CL. By definition, lB positive if the 

P 

principal axis is above the flight path at the nose of the airplane as 
is the case for the curves presented in figure 5 « In general* for 
positive C n p, the expression of R^_ = 0 which does not include 

any C 7 terms is positive and* except for one term* is independent 
£ 

of C n . If* therefore, the coefficient of is positive, R^ = 0 

Is in the first quadrant; whereaB if the coefficient of Is 

negative, R 7 = 0 Is in the second quadrant. As CL increases in a 
positive direction and becomes more positive* the coefficient of 

becomes more negative and R 1 = 0 in the second quadrant shifts upward 
in a destabilizing direction. If Ag Is negative but the absolute value 
of Ap increases* as in going from CL =0.15 to CL = -0.10 in figure 5, 

the coefficient of C, becomes more positive and R-. in the first 

quadrant also shifts upward in a _ destabilizing direction. Thus the results 
indicate that increasing the absolute value of Ag has a destabilizing 
effect on the neutral-oscillatory— stability boundary. 
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According to a previous discussion herein, variations in 

maintain A 2 constant cause no shift in the D = 0 "boundary? 
therefore, 13 = 0 is approximated "by D = 0, such c hang es in 


CL that 
When, 



and should have a negligible effect , on the B = 0 boundary. In 

order to .test this point, calculations were made for a free-flight airplane 
model for Cn varying from O. 3 O to O .63 while s imul taneously 

p *• 

varying in order to maintain the same positive value of A 2 . The 

results showed the expected insensitivity of the E = 0 boundary to 
these changes. 


It should be remembered that D = 0 in the first quadrant is the 
neutral-oscillatory— stability boundary for the long— period oscillation; 
and if instability were to occur, the pilot might not find this type 
of instability difficult to control. 


Effect of C nr , 



and Kj on the Branch of E = 0 


Approx im ated by D = 0 


The D expression indicates that the D = 0 boundary is inde- 
pendent of the derivatives C n ^ and Cy^ and the mass parameter Ey. 

Figures 6 and 7 show a comparison of the results obtained by the complete 
calculations with D = 0 for the cases in which and Cyp, respec- 


tively, were arbitrarily doubled in value. As noted in the figures, Cj^ 

and Cy^ have a negligible effect on the boundary. The effect of % 

on the branch of E = 0 which may be approximated by D = 0 is shown in 
figure 8 . Complete calculations -were made to obtain the E = 0 curves 
for the previously discussed free— flight airplane model. The value of Ag 


used in these calculations was 0.17* For purposes of comparison, E-y was 


arbitrarily increased by a factor of 2.5* Again the results show 
practically no effect of Ej on this branch of E = 0, as is indicated 

by the simplified expression D = 0. For the case discussed in figures 6 
to 8 , the branch of E = 0 approx i mated by E]_ = 0 is in the second 


quadrant and has little practical importance. Hence, the effect of these 
parameters on Ej_ waB not determined for this particular case. 
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Effect of Product of Inertia on the Branches of R = 0 
Approximated by Rj_ = 0 and E = 0 

The product of inertia has been shown to have a "very pronounced effect 
on the lateral stability of present-day airplanes designed for high-speed 
high-altitude flight (references 1 and 7) • The importance of the product 
of inertia is illustrated in figure 9(a), which presents the R = 0 
boundaries of the hypothetical airplane discussed in reference 1 for two 
angles of inclination of the principal axi s relative to the flight 
path, T) = 0° and tj = 2°. Calculations were also made for these cases 
using R = 0; and the results presented in figure 9(a) show the same 

marked stabilizing shift in the boundary, caused by the 2° inclination 
of the principal axis above the flight path, as obtained by the complete 
calculations. The value of A^ for the R^ = 0 calculations was -0.l8. 

The value of c y,^ was then increased so that Aq was equal to 0.13 

(fig. 9(b)) . In this case, D = 0 appears in the first quadrant 

and R 1 = 0 is in the second quadrant. Although both D = 0 and R^ = 0 

are valid boundaries, the results are discussed only for the effect of 
product of inertia on D = 0 since only the combinations in the 

first quadrant are usually of practical significance. Calculations 
for T) = 0° and q = 2° were made using D = 0 and R = 0. Although 
the product-of— inertia factor gyg does appear in the D expression 

(in the term -2 Crffsiy.), 811 examination of D indicates that this term 
could have only a negligible effect on I) = 0 when C^ is much greater • 

than 2C T K y7 , as is usually the case. Eigure 9(h) shows that the results 

predicted from D = 0 agree very well with the results obtained from the 
complete calculations. 

Effect of Radii of Gyration on the Branch of R = 0 
Approximated by R^ = 0 

Figures 10 to 12 are presented for the purpose of showing the close 
agreement between results obtained by using R-j_ = 0 and results obtained 

from reference 1. The three figures illustrate the effect of the radii 
of gyration in roll and yaw kj^ and kg q , respectively, on the neutral- 

oscillatory— stability boundary. Figure 12 emphasizes the fact that the 
simplified expression is sufficiently accurate to predict the effect 
of k T on the oscillatory— stability boundary throughout the entire 

-&o 

range of variation of k^. 



MCA TN No. 1727 


23 


Effect of Wing Loading and Altitude on the Branches of B = 0 
Approximated by Bj_ = 0 and D = 0 

The effects of wing loading and altitude on the neutral-oscillatory- 
stability boundaries were determined simultaneously by consider ing 
variations in the relative density factor because varies 

directly with both wing loading and altitude. - 'An examination of the 
expressions E^ = 0 and D = 0 Indicated that increasing causes 

a slight destabilizing shift in Ej_ = 0 but does not affect D = 0 
since does not appear in the expression for D = 0. The trend Bhown 

by these results agrees with the results found in reference 1 concerning 
the effect of p t on the neutral-oscillatory-stability boundary. 


Comparison between Neutral-Oscillatory-Stability Boundaries 
Obtained by Exact and Simplified Expressions for a 
High-Speed Experimental Airplane 


Some of the neutral-oscillatory— stability boundaries obtained from 
recent calculations for several experimental high-speed airplanes have 
appeared much different from the conventional stability boundaries. 

Because of the complexity of the complete expression for B = 0, it 
is difficult to determine the reasons for such unusual looking curves 
and the significance of the boundaries. From the simplified expressions, 
however, a complete analysis of the boundaries can be easily obtained. 

The B = 0 boundaries of an experimental airplane are shown in 
figure 13(a). In addition to the E = 0 boundaries, the D = 0 boundaries 
are also plotted in the figure. As mentioned at the outset of the 
paper, B = 0 is a neutral-oscillatory— stability boundary only if D 
is positive. The . E = 0 boundaries on the hatched side of D = 0 are 
not therefore neutral-oscillatory— stability boundaries. The boundaries 
for the same experimental airplane calculated from the simplified expres- 
sions are plotted in figure 13("b). The Bj_ = 0 and D = 0 boundaries ~ 

which are not neutral-oscillatory-stability boundaries, as determined by 
the analysis presented in the section entitled "Validity of D = 0 
and B 1 = 0 as IT eutr al-0 s c i llat ory— S t ab i li ty Boundaries," are shown as 


dash-line curves in the figure. In D = 0, the coefficient of C 7 

becomes zero at = O.O56 and, therefore, the D = 0 boundary 

approaches positive infinity in the second quadrant at C n = 0.05 6. 

P 

As C n increases above O.O56, D = 0 returns from negative infinity 
p • 

and appears in the first quadrant. Similarly, B^ = 0 approaches 


negative infinity when CL, is approximately equal to 0.25 since the 

p 

coefficient of in B^ = 0 ^ A^ K yi/ — A^A.^ zero a "^ ‘this value 

of C n . Above C n of 0.25 j B^ = 0 returns from positive infinity 

P p 


Above 
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and appears in the second quadrant. It is necessary to note that in 
figure 13 (a) the neutral-oscillatory— stability boundary is one continuous 
curve; whereas in figure 13 (b) this boundary is composed of two sections, 
one section of Eq = 0 and the other section of D = 0. The latter fact 

provides the important information that the period of the oscillation 
which becomes unstable upon passing through the D = 0 boundary is 
comparatively longer than the period of the. oscillation which becomes 
unstable upon passing through the Bq = 0 boundary. 


CONCLUSIONS 


The following conclusions were reached from a theoretical investi- 
gation of a simplified method for obtaining and analyzing the neutral- 
lateral— oscillatory— stability boundary: 

1. A necessary condition for the lateral-neutral— oscillatory- 
stability boundard is that E = D(BC - AD) - B^E = 0, where A, B, C, D, 
and E are the coefficients of the lateral-stability equation. The 
expression for E = 0 is approximated by the expressions D = 0 

and Bq = BC — AD = 0. Criterions are derived which, if satisfied, 

indicate that the approximate expressions satisfy the necessary and 
sufficient conditions for a neutral-oscillatory— stability boundary. 

2. If D = 0 and Eq = 0 approximate B = 0, the curve D = 0 

represents the neutral-oscillatory— stability boundary for the oscil- 
lation which has a period comparatively longer, than the period of the 
oscillation for which Bq = 0 is the boundary. 


3 . In general, the results of the computations obtained from Bq = 0 

and D = 0 show very good agreement with the results calculated by the 
exact expression for E = 0. Specifically, the results of the 
investigation indicated: 

1 

(a) An increase in the absolute value of the parameter Ag, 
which is equal to 

coefficient due to rolling- : angular— velocity factor, is the trim 

lift coefficient, and Kg is the radius of gyration in yaw ^ causes 
a destabilizing shift in the branches of E = 0 approximated by D = 0 
and Eq = 0. 


2 C lV 


("where CL, is the yawing-moment 


(b) The branch of E = 0 approximated by D = 0 mainly depends 
upon the parameter Ar, and the damping— in— roll derivative • The 

P 

product— of— inertia term also appears in D, but it has a 

negligible effect on. the branch of E = 0 approximated by I) = 0. 
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(c) An increase in the relative-density factor causes a 

destabilizing shift on the branch of R = 0 approximated by R^ = 
but does not affect the branch of R = 0 approximated by D = 0. 

4. The neutral— oscillatory— stability boundaries computed from, the 
simplified expressions show excellent agreement with the correapondl ng 
boundaries presented- in NACA TN No. 1282. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., August 4, 1948 
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APPENDIX A 

DERIVATION OP TEST FUNCTIONS Ar AND Ad 


For a given value of C. 


selected from tile probable range of C r 


y wvAov uvu. xx v/xu. uxuj ^/x ^ ls umxu x c.m if«,o vx Vm 

P 0 

of the airplane for which the lateral— stability ana l ysis is to he made, 
let 


% = r l C 7 p + r 2 


Sr, 


BC- 


= r l 


(%) 


- r r 


= r 


E l=° - 1 - 


D = d i C 7 p + <*2 


SD 


bC- 


= dn 


(%JU 


2 _ 


= d 


B 2 E = ei C z + e 2 


ac. 


= e 


— 6r 


( 2 p)b 2 E=0 ®1 


= e 


As shown in figure lh the exact roots of R = R n D — B^E = 0 occur at 

2 

the intersection of the straight line B E with the parabola R^D. In 
the vicinity of the point C-^ = r, at which R^ = 0, the curve RjD 

is approximated weld by a straight line tangent to the curve at ^7p ~ r '* 
that is. 


RpD 


/ cJR-,d\ . 

~ hy (% - 0 = radl + rid2 X% ~ r ) 

\ Vr 1=0 


If there is a root of R = R^D — B 2 E = 0 near R-^ = 0 (that is, if B 2 E 
intersects R-jD near the point r, in fig. 1^) then 

R = (- r 2 d 1 + - r) - e-^ - e 2 = 0 
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Thus, the approximate deviation of a root of B = 0 from B^ = 0 Is 
given "by 


Ar 


r l e 2 ~ r 2 e l 


/ \ x t. c. a. 

= V ■ l P J - r^dg - r 2 d 1 - e x ) 

9l (e - r) 
r 1 d 1 (d -r) + e ± 


(Al) 


If this deviation, Ar, is small, then B 1 = 0 1 b a close approximation 
to one branch of E = 0. A suitable criterion for this approximation is 


Ar 

< 

r 





5 


or 


Ar 


= 0.01 


whichever is the larger. 

In the case of D = 0, a similar analysis results in the test 
function 


Ad = 


e L ( e ~ d ) 

r l d l( r “ d ) + e l 


(A2) 


If Ad Is small, D = 0 may then be considered a close approximation ’ 
to the other branch of E = 0. A sui' table criterion for this 
approximation 1 b 


Ad 

< 

d 





5 


or 

I Ad I = 0.01 


whichever is the larger. 
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The expressions for r.^ d-^ and e^ for use in equations (Al) 
and (A2) are 

r l = ~ vs) 

d l = ^-*2 


e =2^20^ 


where 


A 1 = *T C tl t + + C 7 p K Z 

^2 = C n p “ 2C L K Z 2 


*3 *^V“^z 2 
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Figure 1.- System of axes and angular relationships in flight. Arrows 
indicate positive direction of angles, t) = e - 7 - e . 
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Figure b.- Fffect of and A2 on the neuirai-osciiiatory-smbiiity boundary 






Figure 6.- Effect of C n on the neutral -oscillatory -stability boundary. 











Figure 8.- Effect of K x on the neutral -oscillatory-stability boundary. 






(a) Negative A£. 

Figure 9.- Effect of K xz on the neutral -oscillatory-stability boundary. 






(b) Positive A2- 
Figure 9.- Concluded. 
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(a) R = R 1 D - B 2 E - 0. 

Figure 18.- Oscillatory stability boundaries for an experimental high-speed airplane. 
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